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We review the anomalies of the HMF model and discuss the robusteness of the glassy 
features vs the initial conditions. Connections to Tsallis statistics are also addressed. 

§1. Introduction 

In this work we discuss glassy dynamics and nonextensive thermodynamics in 
the context of the so called Hamiltonian Mean Field (HMF) model, a system of in- 
ertial spins with long-range interaction. This model and its generalizations has been 
intensively investigated in the last years for its anomalous dynamical behavior. ^^^^^^ 
With respect to systems with short-range interactions, the dynamics and the ther- 
modynamics of many-body systems of particles interacting with long-range forces, 
as the HMF, are particularly rich and interesting. 

The Hamiltonian Mean Field model is exactly solvable at equilibrium, but ex- 
hibits a series of anomalies in the dynamics, as the presence of quasistationary states 
(QSS) characterized by: anomalous diffusion, vanishing Lyapunov exponents, non- 
gaussian velocity distributions, aging and fractal-like phase space structure. Thus, 
it represents a very useful "laboratory" for exploring metastability and dynamical 
anamalies in systems with long-range interactions. The model can be considered 
as a pedagogical model for a large class of complex systems, among which one can 
surely include self-gravitating"'^'^) and glassy systems, but also systems apparently 
belonging to different fields as biology or sociology. Similar features were recently 
found also in the context of the Kuramoto model, ^'^^ one of the simplest models for 
synchronization in biological systems.^''-* Moreover, the proliferation of metastable 
states in the vicinity of a critical point in the phase diagram seems to be quite a 
general feature . 

In this paper we focus on two different aspects of the HMF model: its glassy- 
dynamics and the possible connections with the generalized thermodynamics. We 
will study in particular the dependence of these features on the initial conditions 
considered. 
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§2. Out-of-equilibrium dynamical anomalies 

The Hamiltonian of the HMF model is given by 

N 2 1 ^ 

H = K + v = Y,^ + ^Y.^^-cos{e,-ej)] . (2-1) 

i=l i,i=l 

This Hamiltonian describes classical XY-spins (or rotators) with unitary mass and 
infinite range coupling, but it can also represent particles moving on the unit circle. 
In the latter case the coordinate 9i of the i — th particle is its position on the circle 
and Pi is its conjugate momentum (or angular velocity). The potential energy term 
is rcscaled by in order to get a finite energy density in the thermodynamic limit 
N ^ oo Associating to each particle the spin vector ~Si = (cos ^j, sin^j) , one 
can introduce the following mean-field order parameter M = jj \ X^i^i j which 
is the modulus of the total magnetization. 

At equilibrium the exact solution predicts a second-order phase transition. The 
system passes as a function of energy from a low-energy condensed (ferromagnetic) 
phase with magnetization M ^ 0, to a high-energy one (paramagnetic), where the 
spins are homogeneously oriented on the unit circle and M = 0. The caloric curve, 
i.e. the dependence of the energy density U = H/N on the temperature T, is given 
hy U = ^ + ^ (l — M^) . The critical point is found at a temperature Tg = 0.5, 

which corresponds to the critical energy density Uc = 0.75. 

At variance with the equilibrium scenario, the out-of-equilibrium dynamics shows, 
just below the phase transition, several anomalies before equilibration. More pre- 
cisely, starting from water bag initial conditions in the momenta, i.e. velocities 
uniformly distributed, and a fully magnetized state, i.e. all the = so that 
M(0)=1, the results of the simulations, in a special region of energy values (in par- 
ticular for 0.68 < U < Uc) show a disagreement with the equilibrium prediction for 
a transient regime whose length depends on the system size N. In such a regime the 
system remains trapped in metastable quasi-stationary states (QSS) with vanishing 
magnetization at a temperature lower then the equilibrium one and shows strong 
memory effects, correlations and aging. Then it slowly relaxes towards Boltzmann- 
Gibbs (BG) equilibrium. This transient QSS regime becomes stable if one takes the 
infinite size limit before the infinite time limit^^ . 

2.1. Hierarchical structures and glassy dynamics 

In order to investigate the nature of this trapping phenomenon, it is interesting 
to explore directly the microscopic evolution of the QSS. A way to visualize it is by 
plotting the time evolution of the Boltzmann fi-space: each particle of the system 
is represented by a point in a plane characterized by the conjugate variables 9i and 
Pi. It has been shown^) that, during the QSS regime, correlations, structures and 
clusters formation in the ^u-spacc appear for the M(0) = 1 initial conditions (Ml ic), 
but not for initial conditions with zero initial magnetization (MO ic). In fact, while 
the former case is characterized by a sudden quenching from an high temperature 
state that causes a strong mixing of particles in the metastable regime, in the latter 
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Fig. 1. For the energy density U — 0.69, A*' — 10000 and for different initial conditions we plot a 
snapshot of the jj, — space at time t = 200 inside the plateau regime. The figure illustrates four 
different initial magnetization cases, M=l (top-left), M=0.95 (top-right), M=0.8 (bottom-left) 
and M=0 (bottom right). Competing clusters of different sizes are clearly visible only for the 
case M=l. See text for further details. 



case no quenching is present: both the angles and velocities distributions remain 
homogeneous from the beginning and a very slow mixing of the particles can be 
observed. 

For the Mlic case, the dynamics in /x-space can be also clarified through the 
concept of "dynamical frustration": the clusters appearing and disappearing on the 
unit circle compete one with each other in trapping more and more particles, thus 
generating a dynamically frustrated situation typical of glassy systems. In Fig^ 
we show the dependence of the competing clusters phenomenon as a function of the 



4 



A. Pluchino and A. Rapisarda 



10 



10" 



3 



10" 



:g 10" 

E 
O 



10 



10" 



N= 10000 U=d.69 (a) 

Sums over 20 realizations 



time^200 




10 100 
Clusters Size 



n 



10" 



o 



10 



10" 



s 



10' 



10'' 



• 


M=1 ic 


(b) \ 


o 


M=0.98 ic 




■ 


M=0.95ic 




o 


M=0.9 ic 




V 


M=0.8 ic 






M=0.4 ic 




* 


M=0 ic 





time^SOO 




Fig. 2. We plot for the case U — 0.69 and N — 10000 the cumulative distribution of the clusters 
size, calculated for time t = 200, panel (a), and t = 500, panel (b), in the QSS region. The 
points refers to different initial conditions: while the initial momenta were always uniformly 
ditributed, the initial magnetizations changed as reported in the legend. A power law can be 
observed only for initial magnetization M=l, black circles. A decaying power law with exponent 
— 1.6 is also reported for comparison as dashed curve. See text. 



initial magnetization. The distribution function f{9,p,t) is considered for different 
out of equilibrium initial conditions (i.e.) inside the QSS regime. In particular for 
N=10000 and U=0.69 we plot a snapshot of the /i-space configuration at t = 200 
for four initial magnetizations, namely M=l, M=0.95, M=0.8 and M=0, obtained 
by spreading the initial angles distribution over a wider and wider portion of the 
unit circle. In this way we fix the initial potential energy V{9) and, in turn, the 
magnetization. Then, we assign the remaining part of the total energy as kinetic 
energy by using a water bag uniform distribution for the velocities. 

We define each cluster as composed by particles with both angles and velocities 
in the same /i-space cell. A total of 100x100 cells for the //-space lattice has been 
considered. In Fig^one clearly observes the presence of competing clusters only for 
M=l i.e.. At variance already for M=0.95 i.e. clusters are not very evident and the 
configuration tends to become more and more uniform decreasing the initial mag- 
netization. These simulations show that, although dynamical structures are present 
in the QSS regime up to M=0.4,^'^) frustration due to competing clusters seems to 
depend much more strictly on the initial quenching which is very violent for M=l 
i.e.. In order to give further support to this result, in Fig. I^Jwe show the behavior of 
the cluster size cumulative distributions calculated in the case N=10000 and U=0.69 
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for several initial magnetizations in the QSS regime at time t=200, panel (a), and 
t=500, panel (b). For each one of the 100x100 /^-space cells a sum over 20 different 
realizations of the dynamics (events) has been performed. Then, for each cluster size 
(greater than 5 particles) the sum of all the clusters with equal of bigger size has 
been calculated and plotted. 

As one can see from Fig|2]a power law distribution, which indicates the presence of a 
hierarchical structure, can be observed only for initial conditions with M = 1. In all 
the other cases studied, i.e. for initial magnetization smaller or equal to M = 0.95 
the power law disappears. By comparing panel (a) with panel (b) one can also no- 
tice that, as expected, the distributions does not change significatively in time in the 
plateaux region. In the figure, we report also a power law fit (drawn as a straight 
dashed line above the data points) for the M=l case: for both cases at t=200 and 
at t=500, the cluster size distribution has an exponent decay approximately equal 
to -1.6*). 

Such a strong dependence of the clusters hierarchical structure on the initial con- 
ditions confirms that, although a metastable macroscopic regime is observed in all 
cases, the microscopic structure of correlations is deeply different. The cluster size 
distribution observed for M=l i.e. reminds closely that of percolation at the critical 
point, where a lenght scale, or time scale, diverges leaving the system in a self- 
similar state. ^■^^ More in general, it has been suggested^^-* that, optimizing Tsallis' 
entropy with natural constraints in a regime of long-range correlations, it is possible 
to derive a power-law hierarchical cluster size distribution which can be considered 
as paradigmatic of physical systems where multiscale interactions and geometric 
(fractal) properties play a key role in the relaxation behavior of the system. The 
power-law scaling strongly suggests a non-ergodic topology of the region of phase 
space in which the system remains trapped when the QSS regime is reached. 

These results also supports the so called "weak ergodicity-breaking" scenario 
typical of glassy systems. In ref.^^) such a mechanism has been proposed in order 
to explain the aging phenomenon, i.e. the dependence of the relaxation time on 
the history of the system. ^"^^ Aging has been found also in the HMF model for 
Ml i.e., more precisely in the autocorrelation functions decay for both the angles 
and velocities^) and for velocities only.'') More recently, inspired by the physical 
meaning of the Edwards-Anderson spin-glass order parameter, ^^) we proposed a 
new order parameter for the HMF model, with the aim to measure the degree of 
freezing of the rotators in the QSS regime and to characterize in a quantitative way 
the emerging glassy-like dynamics.^) We define the elementary polarization as the 
temporal average, integrated over an opportune time interval r, of the successive 
positions of each rotator: 

^ f-to+r 

<Ti>= - Ti{t)dt i = l,...,N , (2-2) 

^ J to 

being t^ an initial transient time. Then we average the module of the elementary 

*-* Please nore that due to a misprint, the exponent in the pubhshed version is erroneously 
reported as equal to —1/6. 
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Fig. 3. For the energy density U = 0.69 and diflterent sizes (N=1000, 5000, 10000, 50000), we plot 
the time evolution of the temperature for different initial magnetizations in panel (a) . In panel 
(b), we show the corresponding polarization values obtained for r = 2000 and for increasing 
times along the QSS plateau. The latter seems to be quite constant in time and with a value 
0.24 ± 0.05. We plot a grey zone to illustrate these limits of variation. Notice that for N=1000 
the the temperature start to relax towards the equilibrium value around t=5000, so also the 
polarization tends to increase. 



polarization over the N rotators, to obtain the polarization p: 

1 ^ 

^ = ]^El<^>l • (2-3) 

i=l 

It is important to point out that, in order to calculate the elementary polarization 
of eg 12 -21 one has to subtract the global motion of the system, i.e. the phase of the 
average magnetization, from the phase of each rotator. 

By means of several numerical simulations we showed in ref.^^ that in the QSS 
regime for 0.68 < U < 0.75 and Mlic, the polarization is constant and greater than 
zero - inside the error representing the fluctuations of the elementary polarization 
module - while the magnetization vanishes with the size of the system. Thus the 
polarization can be really considered a new order parameter able characterize the 
glassyness of the QSS regime^^ . 



In Fig|31we present new calculations for U = 0.69 and M = I initial conditions. 
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In panel (a) we report the time evolution of temperature for different sizes of the 
system (A^ = 1000, N = 5000, = 10000 and N = 50000) and after a transient 
to = 1000. In panel (b) we plot the corresponding polarization calculated at different 
times and for a time interval r = 2000. One can see that the polarization remains 
almost constant around a value p = 0.24 and inside an error of ±0.05 (grey area in 
the figure), it does not change substantially inside the QSS temperature plateaux 
and does not depend within the error on the size of the system. In general the length 
of the plateaux regime depends on the size of the system, thus the integration time 
r = 2000 has been chosen in order to make a meaningful calculation for all the sizes, 
included the smaller one {N = 1000). 
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Fig. 4. For the energy density U — 0.69 and N=5000, we plot the time evolution of the temperature 
for different initial conditions for the magnetization in panel (a). In panel (b), we show the 
corresponding polarization increasing the time interval of calculation. The polarization seems 
to depend strongly on the initial conditions and tends rapidly to zero as a function of the initial 
magnetization as reported in the inset. See text for further details. 



In Fig0]we compare the calculation of the polarization obtained for U = 0.69, 
N = 5000 and different initial magnetizations i.e. M = 1,0.95,0.8,0.4,0. We 
report in panel (a) the temperature time evolution, while in panel (b) we plot the 
corresponding polarization as a function of time. Also in this case a time interval r = 
2000 has been considered to calculate the polarization. The plots are in agreement 
with the previous results about the /i-space cluster distributions. In fact also the 
polarization seems to be strongly dependent on the initial conditions and rapidly 
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tends to zero by decreasing the initial magnetization, as reported in the inset of 
panel (b). Again, the fast initial quenching we observe for the M=l i.e. confirms its 
key role in generating glassy behaviour in HMF model. On the other side, as we will 
show in the next subsection, nonextensive features are quite more robust and seems 
to be less sensible to the initial magnetization. 

2.2. Nonextensive thermodynamics and HMF model 

In previous works it was shown that the majority of the dynamical anomalies of 
the QSS regime, among which velocity correlations in the /x-space, are present not 
only for Ml initial conditions, but also when the initial magnetization M{t = 0) is 
decreased in the range (0, 1].^'^^ 

The velocity correlations can be calculated by using the following autocorrelation 
function^'^) 

1 ^ 

Cii) = j;^Y.P^it)pm , (2-4) 

i=i 

where Pj{t) is the velocity of the j-th particle at the time t. In Fig|SJ we plot 
the velocity autocorrelation function for N = 1000, U = 0.69 and M(0) = 

1,0.8,0.6,0.4,0.2,0. An ensemble average over 500 different realizations was per- 
formed. For M(0) > 0.4 the correlation functions are very similar, while the decay 
is faster for M(0) = 0.2 and even more for M(0) = 0. If we fit these relaxation 
functions by means of the Tsallis' q-exponential function 

eg{z) = [l + {l-q)z]T^ , (2-5) 

with z = and where r is a characteristic time, we can quantitatively discriminate 
between the different initial conditions. In fact we get a q-exponential with q = 1.5 
for M(0) > 0.4, while we get q = 1.2 and q = 1.1 for M(0) = 0.2 and for M(0) = 
respectively. Notice that for q = 1 one recovers the usual exponential decay.^^'^^'^^' 
Thus for M(0) > correlations exhibit a long-range nature and a slow power-law 
decay. This decay is very similar for M(0) > 0.4, but diminishes progressively below 
M(0) = 0.4 to become almost exponential for M(0) = 0. 

Velocity correlations seem to be linked to anomalies in the diffusion. In order to study 
diffusion, one can consider the mean square displacement of phases cT^(t) defined as 

1 ^ 

<''it) = j;^Y.i^jit)-(^jm'=<m)-wf> > (2-6) 

where the symbol < ... > represents the average over all the rotators. The 
quantity typically scales as cr'^{t) ~ . The diffusion is normal when 7 = 1 

(corresponding to the Einstein's law for Brownian motion) and ballistic for 7 = 2 
(corresponding to free particles). For different values of 7 the diffusion is anomalous, 
in particular for 1 < 7 < 2 one has super diffusion. We notice that the quantity cr'^{t) 
can be rewritten by using the velocity correlation function C{t) as 

a'^it)= / dti / dt2 < Pj{t2) Pjih) > =2 dti dt2C{t2) , (2-7) 
JO Jo Jo Jo 
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Fig. 5. For the energy density U — 0.69 and N=1000, we plot the time evolution of the correlation 
function C{t) for different initial magnetization. A q-exponential fit is also reported for each 
curve See text for further details. 



SuperdifFusion has been ah'eady observed in the HMF model for Ml initial con- 
ditions.^^ Rercently we have also checked that, even decreasing the initial magne- 
tization, the system continues to show superdiffusion.^") In particular we obtained 
an ananoalous diffusion exponent that goes progressively from 7 = 1.4 — 1.5 for 
0.4 < M(0) < 1 to 7 = 1 for MO .^'^^ No sensitive dependence on the size of the 
system has been orserved. The slow decay and the superdiffusive behavior can be 
connected by means of a conjecture based on a theoretical result found in ref.'^^^ by 
Tsallis and Bukman. In fact in that paper the authors show, on general grounds, that 
nonextensive thermostatistics constitutes a theoretical framework within which the 
unification of normal and correlated (driven) anomalous diffusions can be achieved. 
They obtain, for a generic linear force F{x), the physically relevant exact (space, 
time)-dependent solutions of a generalized Fokker-PIanck equation. Following ref.,^^^ 
it is possible to recover this relationship 

7=J^, (2-8) 



between the exponent 7 of anomalous diffusion and the entropic index q. Hence, 
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being the space-time distributions linked to the respective velocity correlations by 
the eq. 1)2- 7|) . one could think to insert in eq. ()2-8() the entropic index q, character- 
izing the correlation decay and the corresponding anomalous diffusion exponent. In 
order to check this conjecture, we have studied the ratio 2/{3~q) exponent 7 

for various initial conditions ranging from M(0)=1 to M(0)=0 and different sizes at 
U = 0.69 . Within an uncertainty of ±0.1, the data^^-* (not reported here for lack 
of space) show that this ratio is always one, thus providing a strong indication in 
favor of this conjecture. Although rigorous results are still missing the numerical 
evidence here discussed indicate that Tsallis statistics and in particular its entropic 
index q seems to be able to give a quantitative characterization of the correlations 
and dynamical anomalies found in the QSS regime of the HMF model. 



§3. Conclusions 

We have briefly reviewed some of the anomalous features oberved in the dynamics 
of the HMF model, a kind of minimal model for the study of complex behavior in 
systems with long-range interactions. We have also discussed how the anomalous 
behavior can be interpreted within the nonextensive thermostatistics introduced by 
Tsallis, and in the framework of the theory glassy systems. These two frameworks 
seems to have several links which will be further explored in the future. 
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